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Abstract 

The generalization of the effective action [1] of the quark-antiquark system in the 
confining vacuum is performed for the case of arbitrary quark masses. The interaction of 
quarks is described by the averaged Wilson loop for which we use the minimal area law 
asymptotics. 

The system is quantized by the path integral method and the quantum Hamiltonian is 
obtained. It contains not only quark degrees of freedom but also the string energy density. 

As well as in the equal masses case [1] two dynamical regimes are found [2]: for 
large orbital excitations (I S> 1) the system is represented as rotating string, which leads 
to asymptotically linear Regge trajectories, while at small I one obtains a potential-like 
relativistic or nonrelativistic regime. 

In the limiting cases of light-light and heavy-light mesons a unified description is 
developed [2]. For the Regge trajectories one obtains nearly straight-line patterns with 
the slope very close to l/2iro and 1 /na correspondingly. The upper bound on the light 
quark(s) masses which doesn't change considerably this property of the trajectories is also 
found. 



1 Introduction 



Recently the new approach to study nonperturbative large distance dynamics of quark- 
antiquark system in the confining vacuum has been developed and the Hamiltonian for the 
case of equal quark masses has been obtained [1]. In the present paper we consider the general 
case of arbitrary quark masses. The effective Hamiltonian of the system is derived (for a short 
report see [2]) and properties of its spectrum are analysed. 

With the help of vacuum correlator formalizm [3] we represent gauge invariant Green func- 
tion of the qq system in a form where all dynamics of the interaction is described by the averaged 
Wilson loop operator. Starting from the QCD Lagrangian and making use of the minimal area 



1 



law asymptotics for the Wilson loop we arrive to the string picture of q — q interaction. The 
appearence of the string at large distances has been discussed for a long time [4,5] (for recent 
calculations of the colour-tube or string dynamics see ref. [5] ) and here the world sheet of the 
string coincides with the surface appearing in the area law of Wilson loop. 

According to the vacuum correlators method [6, 3] it is the minimal surface which enters 
the area law asymptotics. This aproximation provides a good calculation scheme at least for 
the consideration of the leading Regge trajectories. 

As in [3,1] the minimal surface is approximated as the world sheet of the straight line 
connecting positions of quark and antiquark with the same time in the meson rest frame. We 
discuss the dynamical motivations of this ansatz for the area, which leads to a local in time 
string interaction usually postulated in the flux tube models picture [5]. In this approximation 
the string (minimal string) may rotate and oscillate longitudinally (stretching and expanding). 

We simplify our problem by disregarding effects due to the quark spins and additional quark 
pairs creation having in mind to come back to them in subsequent papers and perform path 
integral quantization of the system of quark- antiquark with arbitrary masses connected by the 
minimal string. 

It turns out that our effective action can be represeted in a form postulated in accordance 
with the flux tube picture in ref. [5] without direct connection to the QCD Lagrangian. 

We generalize the procedure for determination of the center of masses coordinate as com- 
pared with that of the equal masses case [1] where the symmetry between quarks simplifies the 
consideration. Due to the string contribution to the kinetic part of the action, we are to exploit 
the condition, that the center of masses coordinate is decoupled from the relative one. 

As well as in the equal masses case we find two dynamical regimes of the "minimal" QCD 
string with quarks of arbitrary masses. They are distinguished by the energy- momentum 
distribution of nonperturbative gluonic fields along the string: for large orbital momenta I ^> n r . 
the main part of the energy (orbital momentum) is contained in the rotating string and the 
resulting spectrum in the leading approximation coincides with that of the pure string. For low 
values of I (I <C n r ) the dynamics is described by a relativistic or nonrelativistic (depending 
on quark masses) potential-like approach with almost inert contribution from the string. The 
transition between these two regimes is very smooth and the relativization of heavy quark(s) 
due to increase of l,n r is briefly discussed. 

For two limiting cases of heavy- light and light-light mesons a unified description is developed. 
The resulting spectra in both cases are very simple: practically linear leading Regge-trajectories 
and nearly parallel to them daughter trajectories, corresponding to radial excitations with the 
slopes 1/27TO" and 1/na correspondingly. We estimate also the upper bound on the light quark(s) 
masses, which don't change considerably the slope of the trajectories. 

We note that this pattern for the trajectories is very close to one obtained by numerical 
quantization [5] of the same system which was performed without qualitative consideration of 
the underlieing regimes of quark-string interaction. 

The modification of heavy quarkonia Hamiltonian due to the string contribution to the 
orbital momentum is considered. The corresponding correction to the energy is estimated and 
we find the conditions under which it is comparable with the first relativistic correction. 

The paper is organized in the following way. In Section 2 we derive the effective action of 
the qq system at large distances and discuss the basic approximations, which enable one to 
reduce initial field theory problem to the local three demensional quantum mechanics of quarks 
connected by the straight line string. In Section 3 we represent this action in a gaussian with 
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respect to the quark coordinates form. The Hamiltonian for the "minimal" QCD string with 
quarks is obtained in Section 4 and in Section 5 we analyse the properties of its spectrum. 

In Apperdix A intermediate steps of the transformation to the gaussian representation 
are considered. Technical details for calculation of the first corrections to zero order pure 
potential and pure string approximations for the dynamics are collected in Appendixes B and 
C correspondingly. 

2 Effective QCD action of the qq system at large 
distances 



We start with the Green function for spinless gg-system in quenched approximation, which in 
the Feynman-Schwinger representation has the following form [3] 

G{xx\ y,y)= ds dse' K ~ K DzDz <W(C) > A , (1) 
Jo Jo 

where W(C) is the usual Wilson loop operator 

W(C) =trP exp [ig f A^dzfx] (2) 

with 

K = m\s + \j' Q %{t)dt , i„ = ^ (3) 

and the analogous for K. 

The closed contour C consists of initial and final pieces [x, x] , [y, y\ entering the boundary 
conditions and paths z(t) , z(t) of the quark and antiquark. 

At large distances R > T g ~ 0.2 — 0.3 fm (where T g is the vacuum gluon correlation length) 
one gets from the cumulant expansion for Wilson loop asymptotics (omitting perimeter-type 
terms always coming from the quark-mass renormalization ) [6,3] 

< W(C) > w exp(-aS) (4) 

where S is the area of the minimal surface inside the contour C . 

This asymptotics is well confirmed by lattice simulations and allows one to take advantage of 
the approximate reduction of gluodynamics to the formation of buozeen string between quarks. 
As a consequence we are left at this step with the action depending only on quark coordinates 
z^, z^ (since the minimal surface is determined by the form of these trajectories). 

The reparametrization from proper times t,t to the Euclidean times z = r , z = f 
performed as in [1] yields 

dt=^, di=^ (5) 

%0 z 

which amounts [1] to the following substitution in eq.(3) for kinetic terms 

K^K'= j*[J* +IMl ( T ){l + ?(T)}], (6) 
J 2 ^i(r) 
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dsDzQ — > Djii 

and the analogous for K. Here we have introduced path integration over the new functions, 
playing, as we will see, the role of dynamical quark masses 

T T ._ 

M T ) = ^(t) , H2{t) = ^zo(t) (7) 

and T = |(x + x - y - y ). 

Let us discuss the approximate reduction, which we apply to the dynamics of zero compo- 
nents zo , zq. Initial path integral representation takes into account all trajectories, including 
those with the backward motion in time dr = 20(7) dj < , where the signs of zq , zq 
are negative. In the Fock space the backward in time pieces of trajectories are responsible for 
the creation of additional gg-pair. Here we neglect this backward motion of q, q trajectories 
since the backtracking in time quark is dragging with itself also the string and this enlarges the 
action due to the formation of foldings on the string world sheet. 

With respect to the space-time picture of the evolution it means that we don't take into 
account breaking of our string into several ones due to the quark pairs creation from the 
vacuum (in addition to the neglect of such breaking due to qq pairs from quark determinant 
which corresponds to the quenched approximation (1) we started from). 

In what follows this no-backtracking time approximation leads to reduction of quark dy- 
namics to that of the valence quark (connected by the string) Fock sector, where the conditions 

/h(t) > , /i 2 (r) > (8) 

are valid and the transformations (5) are unique. 

Now one is to find reasonable ansantz for the minimal surface of Wilson loop (4) in terms 
of z and z. As in [1,2] we introduce standard compact notations 

f = ( r > P} , 9ab(0 = daW^dbW^, a,b = r,(3 (9) 

so that the area S is 

S = J d 2 Zsfdd~g , (10) 
where w^t, (3) are the coordinates of the string world surface, and 

w.. = , w- 



dr ' " dp ' 

As well as in [7, 1] we can use in the valence quark sector (8) the approximation that the 
minimal surface for given paths z^{r), z^ij) is determined by eq.(10) with given by straight 
lines, connecting points z^{t) and z^(r) with the same r, i.e. we exploit the instantaneous 
approximation to the interaction (10). 

= (r ,z) , z„(t) = (r , z) , (11) 

and 

w im (t,P)^z^t)-P + z^(t)(1-P) , O^P^l. (12) 
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We note that this approximation to the surface geometry is inspired by two limiting cases which 
are of special interest below: in the case / = one can make use of the flat dynamics of quarks 
and in the limit / — > oo quarks and the string are moving along typical trajectories of the 
double helycoid type [7], for which the minimal area indeed is formed by the straight-lines. It 
amounts effectively to the elimination of the second time from the problem and corresponds to 
the instantaneous formation of the string in accordance with a position of quarks. 
Combining the results we obtain the total effective action in the form 

T 1 

A = K' + K' + a J dr J d(3^/detg = K' + K' + A str . (13) 
o o 

where the kinetic terms of quarks are 

T 



+ 



m 2 



+ /i 2 (r){l + l 2 (r)} 



M r)—i vj (14) 

We emphasize, that the presence of the terms K' + K' violates the condition that the ends of 
the string move with the velocity of light. 

We stress [1] also, that it is the valence-quark approximation (8) together with the equal 
time straight-line ansatz (12), which enables one to reduce at large interquark disances nonlocal 
four dimensional dynamics (1) to the local in time three dimensional dynamics (13). 

After integration over /ii, (J>2 our action (13) obtained from QCD Lagrangian coincides with 
that of ref.[5] postulated in accordance with the flux tube picture of the interaction. 

But as we will see our way of formulation and solution of the Hamilonian problem is different 
from numerical one used in [5]. The main advantage of our analitical approach is the possibility 
to reveal underlieing dynamical regimes of quark-string interaction discussed in Section 5. 

3 Gaussian representation for the effective action 

A direct procedure of quantization of eq.(13) is difficult due to the square root term and as in 
[1, 2] we use the auxiliary fields approach [8] to get rid of it. 

The string part of action A str , in (13) is equivalent [1] (for the straight line ansatz) to the 
following one quadratic in quark coordinates [8] 

T 1 

A str , = J dr J dfi V - [w 2 + [ajvfw 2 - 2r](ww') + rfw' 2 ] (15) 



where the field v will play the role of the string energy density and the conditions (11) for 
Zh(t), z^t) are implyed. 

The auxiliary fields v(t, (3) > and rj(r, f3) are integrated out together with /^(t) in the full 
path integral representation for G 

G = J DzDzDvDr]Dii l Dn 2 e- A . (16) 
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We emphasize that (up to the preexponential factor, which is immaterial in our case for the 

derivation of effective action) the integration over //i,// 2 , f and rj effectively amounts to the 

replacement of them by their extremum values [1]. 

In the valence quark Fock sector (8), (12) one can introduce relative r M and center of masses 
coordinates in a selfconsistent way and fix the rest frame of the meson. Contrary to the 

equal masses case [1] here one is to use explicitly the condition that R^ is decoupled from r M . 
First we note that the action (13), (15) must not contain an explicit dependence on R^ (for 

the total momentum P M to be conserved). Therefore the relative coordinate is to be defined as 

r » = < = ^(r) - z^t) . (17) 

Second, to diagonalize the quadratic in z , z kinetic part of A in such a way that the 
Jacobian 

- 1 (18) 

we are to introduce 

R, = C(t)z»(t) + (1-C(r))^(r) (19) 

where the parameter £(r) is determined from the condition, that R is decoupled from r. 
Actually the action can be rewritten in terms of i? M ,r M ,r M as 

^ r 2 2 i 
A = ( dr -^1 + ^1+ \ ai R 2 + 2a 2 (Rr) - 2 Cl (Rr)- 

-2c 2 (fr) + a 3 r 2 + a 4 r 2 }] (20) 

where we have used = R^ + (P — C)r M to express w in terms of R, r and the following 
notations have been introduced 

i i 

ai = J dPbn +H2 + V), a 3 = Jd(3 (^(1 - C) 2 + /i 2 C 2 + (13- Cf") , 



a 2 = I dp (in - C(/ii + /i 2 ) + (P~ C» , a 4 = fd[3(^ + rfu) , 

o v 7 

1 i 

ci = / d(3r)v, c 2 = J df3r}((3 - ()w ■ 

o o 

For the diagonalization of the kinetic part of A one is to put a 2 = 0, which leads to 



(21) 



c{ x = »i(t)+ < (3 > J v(t, [3)d[3 

U ] Ln(T)+ii 2 (T) + fv(T,(3)d(3 { ] 

where 

fpvdp , , 

<f3> =W' (23) 

It should be stressed, that eq.(22) for £(r) which determines the -center of masses" coor- 
dinate (19) involves not only quark dynamical masses ^ but also the dynamical string energy 
density v(t,P). This corresponds to the string contribution to the kinetic part of action (20). 
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Eqs.(19), (22) generalize the equal masses case mi = m 2 [1] where one gets for the extremal 
values of fa 

Mr) = Mr) . (24) 
Since the extremal values of the function v{(3) (that only contribute to the action) are even 



under the exchange (3 



— {13 — |) (due to the symmetry under permutation of the ends of 



the string) one obtains from (22), (24) for mi = m 2 . 

J {(3- l/2)ud(3 = , 



COO = 1/2 , 
R(r) = 1 (z(r) + I(r)) . 
In the end of this section let us fix the frame of reference. For this purpose we perform 

— * — * 

usual canonical transformation from R to the total momentum P 



(25) 

(26) 
(27) 



PR- H(P,...)\dr 



J DRexp i J L(R,...)dr = J DRDPexp 



where H(P, ...) = PR — L(R, ...) and conserved total momentum is defined as 

p = dL(k,...) ^ 
OR 



(28) 



(29) 



Meson rest frame correspons to P = 0. Taking into account conditions (11) leading to 

R (t) = 1, r (r) = (30) 
we obtain for the effective action A, equal to H(P = 0, . . .) 



1 T 



dr 



1 2 1 

m\ wij 1 

Hi fl2 CLl 



{a 3 aif 2 - 2c 2 ai(rf) + (a 4 ai - c\)r 2 + a? j 



(31) 



4 Derivation of the Hamiltonian for the "minimal" QCD 
string with quarks 

Action (31) contains (through the dependence of a iy q) the auxiliary functions fj,± , fj, 2 , v , V ■ 111 
contrast to the first three ones (playing the roles of quark dynamical masses and string energy 
density respectively as one can see for example from eq. (22) ) the function t)(t ,(3) is an 
intermediate one. To derive the effective Hamiltonian for the minimal QCD-string with quarks 
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we perform gaussian integration (A is quadratic form in if) over 77 to obtain the effective action 
A' 

J Drj(T : P)exp[-A] ~ e~ A ' . (32) 
ft amounts to the substitution into A the extremal value of 77 (see Appendix A for details) 



Vext(r,P) 



[rr) 



P 



Hi 



fii + lii, 



(33) 



which leads to the following definition of A' 



x 



[r x f\ 2 _ {fr) 



_l_ r 



c? 2 



a 2 d/3' 



r a*a 
J v 



(34) 



where 



fi 



111 + A*2 

Expression (34) forms the basis of our further calculations. 

To perform canonical transformation from f to p we separate longitudinal and transverse 
with respect to r components of r 

'r 2 = -L{(*r ) 2 + [fxff} (35) 

and obtain for the logitudinal and transverse components of the momentum respectively 

{pr ) 2 



Vr = 



~2 \ rr )_ 



(36) 



^ 2 [pxf 12 
PT = — =TT 



= L(i-C) 2 + ^C 2 + / d^-C?v 



[fx f] 2 



(37) 



We note that the string energy density v doesn't contribute into the longitudinal component 
of the momentum, because the momentum density of the string 



pstr 



w 2 w' 2 — (ww 



') 2 ) 



is orthogonal to w' = r M 



(>*V) = (P str r) = . 



(38) 



(39) 



The standard derivation of H from the action (13) yields (in the Minkowski space-time) for 
the Hamiltonian in functional integral 



H(p,r, u, Hi,ii 2 ) = 2 



(p 2 + m\) (p 2 + ml) _ 2 

— + — + fii + (I2 + r 2 



Hi 



fl2 



a 2 d(3 



v 
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+ 



1 

J vdf3 



+ 



L /r 2 

Vi(i-C) 2 + /i 2 C 2 + )^(/3-C)V 

v. i 



(40) 



where L = p T 2 f 2 . 

This is the resulting Hamiltonian for the "minimal" QCD string with quarks. For the case of 
equal masses mi = m^, which amounts to the possibility of substitutions (24) - (27), it has been 
derived in ref. [1]. Approximate expressions for eq.(40) with disregard of string contribution to 
the orbital momentum were obtained in [7]. 

Let us stress again, that eq.(34) contains auxiliary fields fj,i, fj,2, v. After integration over 
them in the full path integral representation (16) only saddle points values fi^ xt (r), u ext (T,f3) 



dH 

dfii{r) 



o, 



1,2 



5H 



Mr, P) 







(41) 



contribute [1] in the effective Hamiltonian. The conditions (41) lead to the following equations 
for the extremal values of auxiliary fields 



-» 2 i 2 

p z + m{ 

1WT 



-Q i 2 

p + m.2 



a 



1(1 + 1)/? 



l-/(/+l)/f 2 (^-l N 

\ tt 3 P2/ 



= 1-1(1 + 1)/? 



(42) 

(43) 
(44) 



u 2 (t,(3) u 3 

where a 3 = fii(l - C) 2 + + I df3(f3 - () 2 vd(3 and C is defined by eq.(22). 

Only after substitution of these extremal values into the path integral Hamiltonian (40) one 
is to construct (performing proper Weil ordering [9]) the operator Hamiltonian acting on the 
wave functions. 



5 Dynamical regimes of the "minimal" QCD string with 
quarks 

In this section we perform analysis of dynamics of the "minimal" QCD string with quarks. First 
we consider the case when both (heavy) quarks are nonrelativistic. We calculate the string 
contribution to the orbital momentum of heavy quarkonia and find the domain of quantum 
numbers I, n r , where the corresponding correction to the energy is comparable with the leading 
relativistic correction. Second we analyse another limit of one or both light quarks (m 8 < ^Ja) 
and find two relativistic regimes for large and small orbital momentums. Finally we discuss 
finite mass effects and the transition from nonrelativistic dynamics of heavy quark(s) to the 
relativistic one. 

A. The limit of nonrelativistic potential dynamics 
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Let us consider the heavy quarks limit of the QCD string with quarks 

mi, m 2 > \[cr . (45) 

In order to deal with a realistic case we impose also the condition that the characteristic 
distance < r > between heavy quarks is more than the vacuum gluon correlation length T g ~ 
0.2/m and the dynamics is governed by the confining potential. It amounts to the following 
constraint on the masses 

< r > ~ (2ma)- 1 / 3 > T g , (46) 

where m is the reduced quark mass. 

In this case one easily gets in the leading order from eqs.(42)-(44) 

fif Bt = m i , i = l,2; 
v ext = o\r\ , 



m l + m 2 

with the second term in righthand sides of eqs.(42)-(44) being much less than all others. There- 
fore we obtain neglecting all relativistic and string corrections the nonrelativistic Hamiltonian 
in the form (for large distances between quarks) 

H = m 1 + m 2 + -Vp 2 + a\ r\ , (47) 
2m 

where 

m 



mi + m-i 

The first corrections to eq.(47) can be easily obtained from the eqs.(40)-(44) with the result 
(for equal masses m x = m 2 = m case) 

AH- = -14 - 1 <L±?Z . (48) 
4 m 6 bm z r 

The first term comes from the expansion of the ordinary square root 2y/p' 2 + m 2 and the second 
one corresponds to the string contribution to the orbital momentum (see the next subsection 
for a more detailed discussion). 

Simple estimates (see Appendix B for details) give the following results for AE rel 

where E^ t is the eigenvalue of the nonrelativistic Hamiltonian (47) after subtraction of heavy 
masses. One can represent in the form [7] 

EZ = (2m)-^a nr , l , (50) 

where a nr j is the eigenvalue of the corresponding dimensionless Hamiltonian. Taking for ex- 
ample n r — 0, 1 — 1, a ,i = 3.36 [7] we obtain that for not very large I ^ 0,n r the string 
contribution (the second term in eqs. (48), (49) ) is of order of the first relativistic correction. 
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Making use of the asymptotical expression [7] 



a nr ,i 



n 



2/3 / /Or? 

3 + ^ 

\ n 



I » n r , 



(51) 



where n = n r + / + 1 we conclude that for large I > n r these corrections remain to be comparable 
while for large n r ^> / the relativistic one becomes obviously dominant. 

In conclusion we note that the string correction to the nonrelativistic Hamiltonian (47) was 
also considered in ref.[10]. 

B. Relativistic potential regime of the system for small orbital momenta 

Let us consider relativistic dynamics of our system and prove the existence of two relativistic 
dynamical regimes of the "minimal" QCD string with quarks: potential one for small orbital 
momenta (or n r ^> I) and the string one for large I ^> n r , which join each other very smothly. 
For the equal masses case mi = rri2 this analysis has been performed in ref.[l] and here we 
generalize it for the system with arbitrary masses of quarks. 

We begin our consideration of the general case with derivation of the relativistic Hamiltonian 
for zero orbital momentum 1 = 0. As we shall prove it provides a good zero approximation for 
not very large / and for excitations with radial quantum numbers n r ^> I. 

For the case / = eqs.(42)-(44) give 



Hi{r) = \jp 2 + m 2 , 



1,2, 



u{t, ft) — a | f 
so that the Hamiltonian (40) finally has the form 



H(p, f ) = \Jp 2 + to? + \Jp 2 + ml + a 



(52) 
(53) 

(54) 



where p 2 = p 2 = ^S-j . 

We note again, that in the case of / = the string doesn't contribute into the kinetic terms 
and is responsible for the inert potential term a | f | . 

Expression (54) (valid strictly speaking only for I — 0) is widely used in the context of 
the so-called "relativistic quark models" [11] for arbitrary I. The approximate version of this 
Hamiltonian was derived in ref.[7] . 

Now we consider the case of small values of / or I n r . To develop a perturbation expansion 
for the spectrum it appears convenient to recover in expression (40) the dependence on the total 
momentum p 2 = (p 2 + L 2 /f 2 ) with the result for H 



H(p, r, v, ^ 



p 2 + m\ p 2 + m 2 2 
Aii + A*2 H 1 V 



a 2 df3^ 2 
r + 



Ail 



Ai2 



r o*a 
J v 



+ J d(3v+ { 
o 



L 2 /r 



Ai 1 (l-C) 2 + Ai2C 2 + /^(/3-C)V 
o 



where i = — + — . 



(55) 
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As well as in [1] we will demonstrate now that for not very large values of I (or for I <C n r ) 
the term in curly brakets corresponding to the string contribution to the (orbital) kinetic part 
of the expression (40) 



ff « =L /f 



\ 



Mi-C) 2 + /i 2 C 2 + /^03-C)V * 



can be treated as a perturbation to = H — 
In zero approximation one again obtains 



Hi 



p 2 + m{ , 



a r 



p 2 + m 2 , , 



and 



= \jp 2 + m 2 + \jp 2 + m\ + o | f | 

where p 2 = (p r 2 + L 2 /r 2 ). 

This potential like regime (small /) corresponds (for m 2 < a) to 



(56) 



(57) 



(58) 



< Hi > rsj - < V > 



(59) 



with nearly pure inert contribution from the string. It is important that even for = the 
quark-string interaction generates dynamical quark masses Hi ~ a/o". 

The first order correction is determined by the average of ff ^ over wave functions of 
Hamiltonian (58) and can be represented in the form 



e (i) =< lo*l(l + l) 



where we have introduced 



2 z/(°) 2 



/4W 0) ' 



/2'(°) - V /2(°) 



> 



(60) 



fl- 



ip) 



(61) 



with ^W^^ + I^o). 

To separate problems we first consider two limiting cases: heavy-light meson (hi) with 



Hi = m 2 — > oo 



mi = 



(62) 



and light-light meson (Z.Z) with 

m! = m 2 = . (63) 

where a unifyed description can be developed. Effects connected to the presence of finite masses 
will be discussed in subsection D. 

In these two limiting cases zero order Hamiltonian (58) takes a simple form [11] 



if (°) -m 2 = KJf 2 + a\f 



(64) 
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where K = 2 for /./ mesons and K — 1 for h.l ones (with the subtraction of the heavy mass for 
h.l systems). The corresponding spectra can be expressed in the form [11] 



e$ = , (65) 

where 

{E^) 2 = ^^n r + ^l+ 3 - + 5(n r ,l)"j , (66) 

with X(n r ) = 4 and 5(n r ,l) is a small correction. Neglecting deflection (56) one would get for 
the slope 

5iP = W < 67 > 

and the results of this approximation (for the leading trajectory) are shown as crosses in Fig.l. 
As we shall see deflection (56) from the "relativistic" Hamiltonian (58) gives rise to a change 
~ 25% in slopes of Regge trajectories for (heavy-)light mesons, so that they are very close to 
the asymptotical string slopes 1/Kncr from the very begining of the trajectories. 

To calculate matrix elements of if ^ = H — over the wave functions of Hamiltonian 
(64) we will neglect as well as in [1,2] the dispersion of | p | and | f |, i.e. will perform the 
substitution 

<|p*|™.|f|"> — ► <\p\> m ■ <\f\> n (68) 

In Appendix B it is shown that within the accuracy ~ 10% the value of can be calculated 
in the form 

e ( n! = VKE$> , (69) 

where 

- aH{l + 1} (70) 

^nl ~ (0)3 • U U J 

E n i 

Therefore in this approximation one gets for not very large / 

(e nl -m 2 ) 2 = K(E^+E^) 2 (71) 

(with disregard of the term ~ (e^ ) 2 , since it exceeds the accuracy of the first order calculations) 
and the corresponding predictions for (e n i —m^lKno are shown by open circles in Fig.l. For 
the leading trajectory and not very large / (I < 5) 

e (i)/ e (o) < o.05 , (72) 

with the decreasing effect for the dauther trajectories. But we stress again that as for the 
change in the slope this correction leads to a deflection of the order ~ 20 — 25% as compared 
with pure potential result (67). 

To observe the transition into another regime we consider the dependence of solution (44) 
for z/(r, ff) on p. Taking for simplicity the case m\ = m 2 < \fa (and therefore [i\ = fx 2 , C = 1/2) 
one can estimate the term ~ L 2 in eq. (44) as one giving a correction to of eqs. (57) (in the 
same way as we have done for eq. (60) ) with the result 
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In the limit I ^> n r comparing the unity in right hand side of eq.(44) with the asymptotical 
value of expression (73) 

203-1/2) 2 

we encounter a large deflection (due to increase of the string contribution to the orbital mo- 
mentum) from zero order equation leading to expr.(57). This transition is governed obviously 
by the parameter 

l —r ( 74 ) 

2n r + 3/2 v ; 

To take it into account we are to start for large orbital momenta I » n r with another zero 
approximation for Hamiltonian (40). 

C. Relativistic string regime for large orbital momenta 

Let us consider the case of 



I 3> n r 



(75) 



where as we will prove light quarks in the meson carry only a small part of the total energy 
(orbital momentum) 

< v > > < ^ > (76) 

where % — 1, 2 for 1.1. meson (with m 1; m 2 < y/a) and % — 1 (m x ~ y/a, /j, 2 = m 2 — > oo) for h.l. 
meson. For equal masses case m\ = m 2 it was shown in ref. [1]. 
For / ^> n r one can exploit the quasiclassical condition 



< (r- < r >) 2 > < < r > 2 
and expand [1] the r 2 -depending part of Hamiltonian (40) around the extremum 



(77) 



/(/ + l)(a 2 /^/z/)- 1 



1/2 



o 



(78) 



so that in gaussian approximation 

1 



H 



+ 2 



Pr + mj p 2 r +m 2 2 
/ii /i 2 



/ i 1 (i-C) 2 + /i 2 C 2 + /(/3-C) 2 ^ 



+ A*i + A*2 + J vdf3 + 

2 



(79) 



As in the case of small / we first concentrate on the two limiting cases (62), (63). For 
mi = m 2 the symmetry between quarks allows one [1] to restrict the class of auxiliary fields by 



and consequently 



J {(3-l/2)vd(3 = 
C = 1/2 • 



(80) 
(81) 
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For ji 2 — > m 2 — ► oo one obviously gets in the leading order 

C = , n 2 = m 2 . (82) 

Let perform in the limit (75) the expansion of Hamiltonian (79) in powers of \xjv. In zero 
approximation one gets (neglecting radial dynamics, m 8 and ^ of light quarks) 



(tf (0) - m 2 ) = \ J ud(3 + ( 



0*1(1 + 1) If V /2 
/(/? - (tf - l)/2) W/ 



(83) 



where K = 2 for 1.1. case and X = 1 for h.l. one. We emphasize, that in this approximation 
one recovers the pure straight-line string Hamiltonian [12, 1] without radial excitations, which 
will appear only as a correction to expression (83). 

The extremal value of v^\(3) one can find in the same way as it has been done for 1.1. case 
in [12, 1]. 

m ( AKaJl{l + \)\ V2 1 



This energy distribution along the string leads to the following total energy in zero approxima- 



tion 2 



(e<°>) = Kiray/l(l + 1) . (85) 

Expression (84) for the energy density can be easily interpreted as a rotation of homogeneous 
distribution (57) of potential case if one recognises that factor 

1 1 



i_ (jr(/j-*f!) ! 



represents the standard Lorentz-factor, with v(/3) = K{(5— ^f^) playing the role of the velocity 
of the corresponding elementary piece of the string. 
For 1.1. system 

- v (0) = v (l) = 1 (86) 
and for h.l. system one has one half of the string 

v(0) = , v(l) = 1 . 

Let us find the leading correction to zero approximation (83) and recover our starting 
condition < ^ ><^< v > (for mf < Again we first discuss two limiting cases (62), (63) 

for 

mi = m 2 = , 

or 

mi = , m 2 — > oo . 

In Appendix C it is shown, that expanding eq.(79) in /// v up to the first nonvanishing terms 
one obtains (after integration over jj) the following expression for 

x 2 ) , (87) 



H M = 4k (^H^r 3 /j ((pi + cL im \) 2lz + 
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with 

= (^(rf + !>?)) 1/3 . 

where the sum is over i — 1,2 for 1.1. and i = 1 for h.l. mesons, c = 3l2/j — and we have 

7r<r 

introduced a dimensionless variables 

t^lA = ^ = sfKtf/s^ayVH 1 ' 10 , V l = -d 2 x (88) 

and as before X = 2, 1 for 1.1. and h.l. systems correspondingly. For the case of K = 2 this 
expression was obtained in ref.fl]. 

Within the accuracy ~ 5 — 10% the spectrum of H)}> can be calculated (see Appendix C 
for details) in the form 

e W = VK^ (£) \ (n r + 1/2) 4 / 5 l-^ (89) 
and in accordance with our initial assumption (76) one indeed obtains 

At /V + l/2\ 2/5 



1/ v Z 

in the limit (75) we are dealing with. 

The resulting spectrum for M 2 /Kna, where 



') 



< 1 (90) 



M = e ^+e^ (91) 

is shown by open squares in Fig.l and doesn't depend on K (again we are to disregard the 
term ~ (e^) 2 since it exceeds the accuracy of the calculations). As we argue in Appendix C 
the exact value of (e^) exact ~ 0.9e^ and we take it into account in this figure. It follows from 
Figure 1 that low / and high / approximations join very smoothly forming (within the accuracy 
~ 5%) Regge trajectories with the slope 

dM 2 , s 

—— = Kna (92) 
dl 

and intersepts, determined from the Hamiltonian (64) 

M 2 (l = 0) = Kna ^2n r + ^ + S^j . (93) 

The values for M 2 are very close to results found numerically in ref. [5] (black dots in Fig.l). 
A small deflection between them is out of the accuracy of the first order calculations we have 
done. A short discussion of similarities and differences between our approach and that of ref. [5] 
is postponed till the Conclusion. 

It is important to stress, that " minimal" QCD string combines together properties of both 
the string models (string like slope (92) ) and ones of the potential models (existence of radial 
excitations, which lead to the intersept (93) ). We have found both in the 1.1. and in h.l. systems 
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two dynamical regimes distinguished by the contribution of the string into the (orbital) kinetic 
part of the action: in the potential regime for small I or 

« 1 (94) 



2n r + 1 



the string constitutes inert linear potential with i/(r, f3) close to the homogeneous distribution 
(57); in the string regime for 

> 1 (95) 

the main part of the energy (orbital momentum) is carryed by the string with energy density 
u(t, (3) given in the leading order by expression (84). 



D. Mass effects in the minimal QCD 
string with quarks 

To analyse the dynamics for arbitrary quark masses case we first consider 1.1. mesons with 

mi < y/a , m 2 <^/a (96) 

and h.l. systems with 

mi < y/a , m 2 > \fo . (97) 

Our main aim here will be to estimate upper bound on light masses, which don't change 
considerably Kna slope (92) of Regge trajectories for massless light quark(s), discussed in the 
previous section. We also calculate in this case leading mass corrections to the spectrum. 

As in the previous section we concentrate on two limiting regimes (94), (95). 

In the potential regime (94) zero approximation as well as in the previous section is described 
by Hamiltonian (the subtraction of the heavy mass in h.l. system is implyed) 



= yp- 2 + mf + yjp* + m% + a | f | (98) 
with the spectrum [11] in the case of m <C<| p |>~ E^ 

e (0) « £ {0) (l + ^fj , 5 {0) = VKE<® , (99) 

where E^ is given by (66) and the sum is over i = 1.2 for 1.1. or i — 1 for h.l. mesons. 

We calculate the correction =< H — > in the same way as for the massless case 
(70) of the previous section (see Appendix B for details). 

Making use of the virial theorem and neglecting the dispersion of Hi ,v one gets in the 

leading order of ( ^fa^ 




(100) 
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The substitution of (100) in (60) gives for the correction e^> in massive case 

/gw\ = j(i + i)q» (101) 



Therefore we have (up to the term ~ m 2 ^- ) for the mass squared M 2 = (e^ + e^) 2 

* ( °> 2 (l-^^) 2 + 2Em?(l + 3 



M 2 = K 



1 Ml + l)a 2 , 



E(°) 4 



(102) 



Let us estimate quark mass corrections to the slope of the trajectory. The first term of eq. 
(102) gives in accordance with expression (92) 



dM 2 T ^ 
^ K no . 



dl 

For the leading Regge trajectory the deflection of the slope from the massless case can be 
expressed in the form 



A dM_ a rn 



2 



/(0 £ *- L (103) 



Kira 37T 3 * 
where f(l) = -^r^ with /(l) « 4. 

Hence the slope of leading Regge trajectory is not essentialy changed when 

K o 3 

^ m 2 <— (7 - 40(7 (104) 

where K = 2 for 1.1. and Jl" = 1 for h.l. system. For the daughter trajectories the restriction 
for the masses is even weaker. 

We emphasize that restriction (104) is much more moderate as compared with that for the 
deflection of intersept (93) (the second term in the r.h.s. of eq.(99)) to be small, which amounts 
to 

(105) 

i=i 

To complete analysis let us consider the case / ^> n r when radial quark dynamics can be 
considered as a perturbation [1]. 

We start with expression (87) for the Hamiltonian if ^ which determines the corrections to 
zero (pure string) approximation (83)- (85). As it is shown in Appendix C the string asymptotics 
(where < v >^>< Hu g ht >) is valid for 

> 1 . (106) 



E K 2 

7T<T I 



Keeping this condition satisfied we first consider the case of nonrelativistic radial dynamics 
when p 2 <C J2i m f , which is achieved as we shall see if 



'r 

5 



' » fe) ■ (107) 



(n r + |) 6 \J2im 
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In this limit one obtains (see Appendix C for details) the spectrum of in the form 



It is easy to make sure, that condition (106) is equivalent to l /e < f > <C 1. Making com- 
parison between p 2 , and S^mf one recovers that in the regime (106), (107) radial dynamics is 
indeed " nonrelativistic" 

< vl >~ (£^ 2 M 5 ) 1/6 '- 1/6 K + \) « J2 m t ■ (109) 

We note that condition (109) follows also from the requirement, that the first term of 
eq.(108) is much more than the second one. 

To calculate the spectrum in the opposite asymptotics of string regime (106) 

<vl > > m 2 (110) 

and merge smoothly the massless case (89) one is to consider (see Appendix C) the following 
domain of I, n r 

' « feY . (no 



(n r + |) 6 

where the spectrum of Hjp has the form 

'3\'/ 5 5 



«. = ^ (f) ' l(™) 1/2 r 3/1 °(nr + l/2y /5 + (H2) 



|) 2/6 E"»?(^)- 1/a r 1/10 K + i/2)- 2 / 5 ) 



Restriction (110), (111) comes from the condition, that the first term of eq. (112) (which is 
the massless correction (89) of the previous section) is much more, than the second one. Due 
to condition (106) we again have the restriction e^\/ef^ <C 1 to be satisfyed. 

To conclude this section let us discuss briefly the transition from nonrelativistic dynamics of 
heavy quark(s) (45) in the heavy-heavy or heavy-light mesons at small n r , I to the large n r , I 
regime, when both quarks become relativistic. We concentrate on the connection between this 
transition and the transition from the potential regime of subsection A to the string one. As 
we easily conclude from eq. (66) the relativization, j? ~ m 2 , appears for 

nr + l,) „ 1M + -2), (113) 

7T / 7T a 

If excitation of the system corresponds to n r ^> I , then one achieves the potential relativistic 
dynamics, described by the Hamiltonian (58) with the small string corrections (60). In the 
opposite case of orbital relativization / > n r as it is clear from the discussion of eq. (73) there 
is a considerable diflection from pure potential regime (58). Therefore in this case both quarks 
start to become relativistic in the regime intermediate between potential one (58) and string 
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one (83) (valid under the condition (106) ) where the dependence of v on (5 is considerable but 
still differs from the asymptotical one (84). 

6 CONCLUSIONS 

Let us summarize the results. We obtain the generalization of Hamiltonian [1] for the 
spinless quark and antiquark in the confining QCD vacuum for the case of arbitrary quark 
masses. Starting from the QCD Lagrangian we make use the minimal area law asymptotics for 
the averaged Wilson loop which leads to the appearence of the minimal QCD string connecting 
quarks. The string contributes (for I > 0) in the kinetic part of the effective action and therefore 
in order to introduce the center of masses and relative coordinates we are to use explicitly the 
condition that they are decoupled from each other. Introducing the auxiliary fields we represent 
the action in the quadratic form with respect to the quark coordinates and diagonalization of 
the kinetic part of the action leads to the proper definition of the total momentum P. Choosing 
the rest frame of the meson P = we arrive (neglecting pathes of quarks with backtracking in 
time) at the effective Hamiltonian. Additionally to quark coordinates it contains the auxiliary 
fields playing the roles of the quark dynamical masses /^(t) and the string energy density 
u(t, (3). Integration over them amounts to the substitution of their extremal values after which 
one is to construct the operator of Hamiltonian acting on the wave functions. 

The interaction between quarks and the string gives rise to appearence of two different 
dynamical regimes (as well as for the equal masses case [1]), which are distinguished by the string 
contribution into the kinetic part of the Hamiltonian. For the low orbital momenta (/ <C n r ) 
dynamics is described in the leading order by the relativistic linear potential Hamiltonian with 
almost inert contribution from the string constituting potential a \ f \. In the opposite limit 
of / n r the system behaves as the rotating string which carries the main part of the orbital 
momentum and energy. The transition between these regimes is smooth and relativization of 
heavy quark(s) due to the increase of I is also considered. 

We develope the unifyed description of the heavy-light and light-light mesons and prove 
that in these limiting cases Regge trajectories are nearly straight line with the slope close to 
1/Kna (K — 1, 2 for h.l. and 1.1 mesons respectively). The upper bound on the light mass(es) 
which don't change considerably the slope is obtained and the leading mass corrections are 
calculated for both regimes. It appears that the slope is much less sensitive to increase of quark 
masses as compared with the intersepts of the trajectories. 

We note that our results for the spectrum are very close to that obtained by the numerical 
quantization of the same action [5] postulated without derivation from QCD Lagrangian. We 
stress that as compared with [5] our approach gives the qualitative picture of underlieing regimes 
of quark-string interaction. Also it enables one to avoid in the leading order the complications of 
Weil ordering [9] in Hamiltonian (40). In zero approximation Hamiltonian both for small I (I <C 
n r ) (58) and for I ^> n r (83) is the sum of the terms depending either only on the momentums or 
space coordinates which doesn't require any Weil ordering. The need for this ordering appears 
only for the correction terms (56), (87) and calculation of them can be performed with the 
neglect of the dispersion (Weil ordering). Such approximate procedure gives the value for the 
leading corrections within more than 10% of accuracy. In [5] the numerical procedure makes 
it difficult to keep track of the ordering and relies on the smallness of the dispersion (see 
forthcoming paper for a more detailed discussion). 
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We also consider the modification of the heavy quarkonia Hamiltonian arising from the 
string contribution (at distances > T g ) to the orbital momentum. We estimate the energy 
correction due to this effect and find conditions under which it is comparable with the first 
relativistic correction. 

The authors are grateful for usefull discussions to Yu.S.Kalashnikova, A.B.Kaidalov, 
Yu.A.Simonov, J.A.Tjon and M.G.Olsson. They would like also to thank N.A.Aksenova for 
typing this manuscript. This work is supported by Russian Fundamental Research Foundation, 
grant N 93-02-14937. 
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Appendix A 

In this Appendix we perform the gaussian integration (32) over r/(r, (5) to go over from 
effective action (31) 



i T 



dr 



mi m 



1 



1 1 {03^1^ - 2c 2 ai{rr) + (a 4 ai - cjr +0x1 

/il H2 01 



(A.l) 



to eq.(34) expressed in terms of physical quantities: dynamical masses /ii,/i 2 energy density 
along the string u(t, (3) and relative coordinates. 

In what follows we make use of the fact that gaussian integration effectively amounts to 
substitution of the extremum value of r] ext into eq.(A.l). On the basis of the definitions (21) of 
a,i and q we obtain from (A.l) the following extremum condition for 77. 



2 / d/3u r 2 {a x r] - a) - (rr)on(/3 - C) Srj = 



(A.2) 



that results in 



7/ 



ext 



ci (rr) 
ai r z 



with c± — J df3rj v. 


Resolution of (A. 3) finally gives 



1 ~2 



P 



Expression (A. 4) leads to the following cf x * 



where we make use of 



^ext 



rr 



(C 



/il 



/il + H2 



) Qi 



c 2 



(rr) 



1 

y = ai - (/^i + /i 2 ) , 


1 

J d(3v(3 = (a\ — /ii , 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



_L 

J df3vf3 2 = a 3 - /ii + C 2 °i • 



To calculate (A.l) with rj = rf xt we note that for quadratic in 77 form one has the following 
relation at the extremum (A.2) 



d e 4 xt ai - (cf *) 2 ) r 2 = (rf)df ai , 



(A.7) 
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where we have separated rj- dependent part of a 4 introducing 64 = J rfvdfi. 



With the help of (A. 7) one concludes that ^-dependent part of the action (A.l) amounts to 
the following contribution into A 

- (rr)c e 2 xt a x , (A.8) 
that results finally in the effective action (34) for the rest frame 



1 T r 2 2 1 

A'i^^u) = - [ dr ^ + ^ + / i 1+/l2 + ( ft (l-(f + /l2 ( 2 + /#-()V)x 



[fx f] 2 (fr) 2 „ 
x ^ + ^hpi + r 2 



2 

dp— 
v 



d(3v 



(A.9) 



where /} 



M1+M2 



is the reduced dynamical mass. 
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Appendix B 



In this Appendix we calculate first order correction (56) to the relativistic Hamiltonian (58). 
Making use of the approximation (68), proposed in [1], we neglect in (56) the despersion of 
| p | and | f | (and disregard also Weil ordering [9]) and express =< > in terms 
of < | p | >, < | r | > averaged over wave functions of zero order Hamiltoian (58) and 
connected to of eq.(66). For this purpose we first evaluate < fif^ >, < > (determined 
by eqs.(57) ) for m 2 < y/a and then substitute them to eq.(60) for e^K Also we will develope 
the procedure to estimate the accuracy of this approximation and find that account of nonzero 
dispersion (Weil ordering) would change the result for less then for 10%. 

First we use the virial theorem for zero order Hamiltonian (58) 

( , P2 + i P2 )=*< \ r\ > (B.l) 
\\Jp 2 + m 2 \jp 2 + m\l 

where the average is performed over the wave functions of (58). 

Neglecting the dispersion one easily obtains the system of equations for < ffl >, 
< /4 0) > and < > 

< (0) ml (0) m% (o) 



F 



W =< ^S 0) > + < ^ > + < i/<°> > ( B -2) 



< fi^ > 2 -m\ =< /4 0) > 2 -ml 

where the first one follows directly from (B.l), the second one is just the definition of the 
eigenvalues of eq.(58) and the last equation ensures the condition that the total momentum of 
the meson is equal to zero. 

To obtain from eqs.(B.2) the values of < ^ > and < > with the account of leading 
mass corrections ~ Y^f=\ m \ /c it is sufficient to consider the reduced system of equations in 
two limiting cases of h.l. (62) and 1.1. (63) mesons 

^ ,,(o ^ L>i=i " l i _ ^ v _ > /r> q\ 




^ 1 £(0)2 

where K = 2 for 1.1. system and K = 1 for h.l. one. The solution of the system (B.3) gives 



e (0) 




o ( Sim? 


2K 




6 [e^ 










\ E(°Y 2 / 



(B.4) 



where e(°) = y/KE^°\ defined in (66). 

Substitution of these values into the eqs.(60), (61) results in the following value for 

<if(i)> =e (i) 

e (o) £(o) 4 [ E(°) 2 1 ' 
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where the summation over % runs from 1 to K. 

In the rest of this appendix we estimate the effect of dispersion for the calculation of 
and find it to be less then 10% even for the ground state of the Hamiltonian (58). For this 
purpose we first represent Hamiltonian (58) (taking for the sake of definitness K — 1) in the 
form containing auxiliary fields ji(r) and u{t) 

\ 1/^2 \ 

(B.6) 

Integration over /x(r) and u{r) amounts effectively [1] (as in the main text) to the insertion of 
their extremal values 




/W = V P 2 + m 2 , u ext = a\r\ (B.7) 

which enables to recover eq.(58) for h.l. system. 

To evaluate the effect of dispersion it is sufficient to restrict the integrations DfiDu by 
that over ji and v independent of time. Actually such approximation gives (for not very large 
quantum numbers l,n r ) deflection in the spectrum (and therefore in dispersion) less then 5 - 
10 % (see [3] for the discussion). As a result within this accuracy we reduce initial Hamiltonian 
(58) to one of the oscillator type. Our strategy is to compare matrix elements < |p| n |r| m > 
with their approximate counterparts <| p\> n <\ f \> m where averaging will be performed over 
our oscillator version of eq.(58). For this purpose one is to obtain [3] first the energy E(/i, v) as 
the function of ji and v. After this we find the extremal values of /i and v from the conditions 

dE_ ^ ^ = Q ^ d_E_ ^ ^ = Q ^ (B g) 

and substitute them into the wave functions in order to calculate the matrix elements under 
consideration. 

Along these lines we have for the spectrum of the Hamiltonian (B.6) 

E M = ^ + ^ + (B.9) 

3 

n + 



JJw V 2 



so that 





N(" + §) 



(B.10) 



and finally in zero order of m one obtains the following spectrum 



E^ 2 =4(n + £)a , (B.ll) 

with n = 2n r + / . We note that expression (B.ll) is reasonably close to the exact one (66) 
(K = 1) and even gives the correct Regge slope (67). 

The wave functions of the approximate version of Hamiltonian (B.6) in the momentum 
space have the standard form 

a nlm (p) = Cp l exp i-^f] Y lm (9, <p)F I + ? a 2 p 2 ) , (B.12) 
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where a = ( a2 jf x T ) and C is the normalization constant. 

To take the symplest example let us estimate the dispersion contribution to the evaluation 
of the first order correction Ae ~ (jy) m the relativistic expansion of the energy ef^ for a 
heavy-light meson (neglecting the string contribution (56) to the orbital momentum). Keeping 
only the first term in the expansion one starts with 



=< y/p 2 + m? + a | f\ > « <^ + a\f\ > + ]- < >= ejj? + Ae (B.13) 



-y 777, 

2 



where is given by eq.(64) with K — 1. One is to compare now two quantities 1/ < a/p~^ > 
and < l/^/jT 2 >. 

For this purpose we consider n r = 0, since in this case the discrepancy (for a given /) 
maximal. For such quantum numbers wave function (B.12) is reduced with the help of Cummer 
formula to the simple form 

a n im(p) = C ■ p l exp (-^f^J Y lm {6, <p) (B.14) 
with C 2 = 7^p|y • The dispersion contribution Ae D is defined in the following way 

A£ - < ^ >=m + ¥ (< jj» > -77^) s AE "° + As " (R15) 

with = 2 < //(°) >= 2 < V / P~ ? > given by (B.ll). The difference Ae D between correspond- 
ing matrix elements can be approximately evaluated with the help of wave functions (B.14) 
with the result 

1 _ 1 £0+1) 

V</^>>-<"<°>> ' r(i+2) c B 16 \ 

1 r(*+i) 



~ v /<c t0) >< uW > ' r ( z +l) 

where are defined in (B.10). 

Making use of Vallis formula for gamma functions (which is valid strictly speaking for I >> 1, 
but still gives a required accuracy for / > 1 ) we obtain for the dispersion contribution 



Ae" = ^ ^L— . (B.17) 

^ 8(l+l)(l + i) 

Expression (B.17) as compaired with the Ae D=0 gives the correction to Ae of order of ~ 10%. 
The dispersion correction in the slope occures even smaller (to the analogy with the ratio of 
quark mass contribution (105) in the energy and in the Regge slope (104), discussed in the 
main text). In conclusion of this Appendix we note that relative contribution of the dispersion 
from higher powers of // and v to the energy is also either of the same order or even smaller. 



26 



Appendix C 



In this Appendix we derive expression (87) for Hamiltonian which describes radial 
dynamics in the limit of large I and calculate the corrections to the spectrum of pure string 
zero approximation (85). We exploit the method, developed in [1] and consider in a unified way 
two limiting cases of the full straight-line string with mi = m 2 — m, fii = fi2 = ^ {K = 2) and 
the half of the string with m-i — m, m 2 — > oo (K = 1). We start here with Hamiltonian (79) 
of the main text. Our strategy will be to expand (in two limiting cases under consideration) 
the following part of (79) 



^ M i(l-C) 2 +M2C 2 +/<i/M/3-C) 2 



1/2 



= 2 



1/2 



+ 



(C.l) 
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and then integrate this approximate expression over fi to end up with Hamiltonian (87). 

After expansion (C.l) the terms in eq.(79) independent of fit are collected into the pure 
string Hamiltonian (83) which gives zero approximation to the problem. To obtain the leading 
radial corrections one is to find the extremal value (84) for from eq.(83) (in the way similar 
to that of [12, 1]) and substitute it into H with the result 



H = (Kna^lil + l)) 1 / 2 + \ 



3K 2 fi 2 



+ 



///2 2(Knal) 1 / 2 



+ 



r ) 



2{K7iaiy/ 2 

where the first term corresponds to the pure string result (85). Introducing new variables 

r — r 



(C.2) 



A* 



Kjj, , p = \f~Kp , 



r 



one obtains the unified radial Hamiltonian = H — ei°^ for h.l. and 1.1. mesons 



K 



(C.3) 



H r = -VK 



+ 



3/i 



,~,2 



2(7ial) 



1/2 



+ 



(na) 2 (f 



TO? 



2(%al) 



1/2 



KH r 



(C.4) 



so that H = ef^ + H^> = y/K(E^ + H^). After the integration over jl, which amounts to 
the substitution of its extremal value 



H = 



-(trx/) j 



1/3 



(C.5) 



we arrive at the effective radial Hamiltonian 



3 4/3 [p' + Y, 



yK --i m i 



2/3 



(na) 2 (f-f ) 2 
2(7ra/)V2 



7TO" 



Z)V6 



+ 



(C.6) 
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Introducing instead of (r — r ) a new dimensionless variable x 

(r - f ) _ v* _ 3 2/5 ■ Z 1 / 10 

)l/2 



X 



A/; 



7TO" 



(C.7) 



one can represent in the following form 



' 3 4/5( 



7TO" 



.1/2 



_ r 3/10 



2/3 



(at 



which is convenient for numerical calculations. 

In order to obtain an approximation to the spectrum of we consider (C.4) for the 
restricted class of functions ft, independent on r. This procedure usually gives the accuracy 
about 5% - 10% for the states with not very large n (see also the discussion in Appendix B) 
[1], [3]. To this end [7] we first evaluate for a given Ji 



rn i/2^/2 ^T72 \™r "I- 2y | -I- 2/i 4 ( TO I) 1/2 



(2(7r<r0 1/2 ) 



(C.9) 



The integration over /i, as well as in the general case, amounts to the substitution of its extremal 



value 



dp, 



in (C.9) . The extremal condition has the form 



(2(7T ( T/) 1 /2) 1 /2 £3/2 



trx/) 1 ^ 



(CIO) 



and is difficult to solve analitically. Therefore we find the solution in two different asymp- 
totics. The first one corresponds to the relativistic radial dynamics 



K 



< Pr > > E 



m. 



(C.11) 



i=i 



that in the leading order gives 



H = 



(nam 1 '* ( 1 



1 2/5 



(C.12) 



and therefore the spectrum (C.9) can be represented in the form 



n r ,l 



n r ,l 



with 



(1,0) _ 5 /3V/5 



<f = (!) 2/5 ^- 1/10 (^ + l) 



4/5 



A -2/5 



(C.13) 



We note that is the radial correction corresponding to the massless quark case (84) 

which was evaluated in [1]. 
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To improve the accuracy we take into account that usually the value of e obtained according 
to the time-independent /i ansatz (C.9), (C.10) is 1, 1 times larger than the exact one (see [3] 
for a comparison). Therefore one is to multiply eqs.(C.13) additionally by a factor 0.9. 

The condition (C.ll) can be reformulated as E^f /E^f <C 1 with the result 

(sMY (c . 14) 

V 7T<7 / t 



Together with the condition, that zero order contribution (80), E^ r = \Jna(l{l + l)) 1 / 2 , is 
much more than E^fi of eq.(C.13) 

l/(n r + 1/2) > 1 (C.15) 

one obtains the domain of quantum numbers where expression for the spectrum (C.13) is valid. 
In the opposite nonrelativistic asymptotics for radial dynamics 

<P 2 r> « XX (C16) 



i=l 



one obtains in the leading order 



1/3 



(C.17) 



which gives rise to 



p (l,0) _ (E,:^ 2 ) 2/3 7-1/6 

_ (3) 1 / 6 1-V3 + n 1 J 

^.i - UJ (L>?) 1/6 l n - + 2j • 

As well as in the previous case the requirements of selfconsistency, E^°i ^> E^^, E^l ^> 
E^j, determine the conditions 

rsi^v 5 » <»- + , 1 / 2 ' 6 (c.i») 

y 7TO" J I 

and 

i » (E to ?/^) ( c - 2 °) 

respectively. 

To summarize we calculated radial corrections in different asymptotics of the string regime 
which is valid under conditions (C.15), (C.20) that can be expressed in the unified way 

» 1. (C.21) 
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Figure Captions 



Fig. 1. Regge trajectories of heavy-light and light-light mesons (with m q = for 
light quarks). Crosses are the results of the approximation of eqs.(64) 
- (66) and open circles are calculated with an account of the leading 
correction, eqs. (69) - (71) (they are not shown for I = 5 since they 
practically coincide with open squares). Open squares are the predictions 
of the large I approximation (91) with an account of the correction (89). 
Black dots are the results of the numerical calculation of ref. [5] 
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This figure "figl-l.png" is available in "png" format from: 
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